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Abstract 

Background: Hypoxia-inducible factor (HIF) is a master regulator that mediates major changes in gene expression 
under hypoxic conditions. Though HIF family has been identified in many organisms, little is known about this 
family in schizothoracine fish. 

Results: Duplicated hif-a {hif-laA, hif-laB, hif-2oA, and hif-2aB) genes were identified in schizothoracine fish. All the 
deduced HIF-a proteins contain the main domains (bHLH-PAS, ODDD, and TAD), also found in humans. Evidence 
suggests a Cyprinidae-specific deletion, specifically, a conserved proline hydroxylation motif LxxLAP, in the NODD 
domain of schizothoracine fish HIF-laA. In addition, a schizothoracine-specific mutation was observed in the CODD 
domain of the specialized and highly specialized schizothoracine fish HIF-laB, which is the proline hydroxylation 
motif mutated into PxxLAP. Standard and stochastic branch-site codon model analysis indicated that only HIF-laB 
has undergone positive selection, which may have led to changes in function. To confirm this hypothesis, HIF-as 
tagged with Myc were transfected into HEK 293 T cells. Each HIF-1 ctB was found to significantly upregulate luciferase 
activity under normoxic and hypoxic conditions, which indicated that the HIF-1 aB protein was more stable than 
other HIF-as. 

Conclusions: All deduced HIF-a proteins of schizothoracine fish contain important domains, like their mammalian 
counterparts, and each HIF-a is shorter than that of human. Our experiments reveal that teleost-specific duplicated hif-a 
genes played different roles under hypoxic conditions, and HIF-1 aB may be the most important regulator in the 
adaptation of schizothoracine fish to the environment of the Tibetan Plateau. 
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Background 

Gene duplication was first proposed by Bridges in 1936 in 
a study of the fruit fly Drosophila melanogaster [1]. It was 
further detailed by Ohno in a book, Evolution by Gene 
Duplication, in 1970 [2]. The ancestor of vertebrates expe- 
rienced two rounds of whole-genome duplication which 
may have played a key role in adaptive evolution of verte- 
brate [2-5]. The third round of whole-genome duplication 
that occurred in teleosts is thought to have driven the di- 
versification of teleosts [6-9]. 

Oxygen is critical to all aerobic life, particularly for an- 
imals living in water, which contains only l/30th the 
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oxygen of air at the same partial pressure [10-13]. So- 
phisticated cellular mechanisms have evolved to allow 
organisms to detect and react to changes in oxygen 
levels. In these mechanisms, oxygen acts as the terminal 
electron acceptor in the mitochondrion of metazoans 
[11,14,15]. When studying the expression of EPO from 
Hep3B cells under hypoxic conditions, Semenza and 
Wang [16] identified a nuclear factor bound to the 
hypoxia-inducible enhancer (HIE) of the erythropoietin 
(EPO) gene at a site required for hypoxia activation of 
transcription. This nuclear factor was named hypoxia- 
inducible factor (HIF). In mammals, HIF mediates major 
changes in gene expression under hypoxic conditions and 
participates in a series of processes including angiogenesis, 
erythropoiesis, glucose and iron transport, glycolysis, and 
cell-cycle control [11,17-21]. HIF is a heterodimeric 
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transcription factor of the bHLH-PAS (basic helix-loop- 
helix-per-ARNT-sim) family. It consists of an erratic alpha 
subunit (such as HIF-la) and a steady beta subunit (such 
as HIF-lp 1 or ARNT) [11,19,22-26]. According to studies 
in Caenorhabditis elegans and D. melanogaster, there is 
only one HIF-a in invertebrates, but in vertebrates, there 
are at least three functional HIF-a isoforms: HIF-la, 
which regulates the acute hypoxic response; HIF-2a, 
which regulates the chronic hypoxic response; and HIF- 
3a, which inhibits the activities of the other two isoforms 
[11,14,19,20,27,28]. HIF-la and HIF-2a both have two 
oxygen-dependent degradation domains (NODDD and 
CODDD) that are located in the central region, and 
two transactivation domains — an inner activation do- 
main (N-TAD) that overlaps with the CODDD and a 
carboxy-terminal activation domain (C-TAD). HIF-3a 
lost the C-TAD [14,19,22,29]. With adequate oxygen 
(normoxia), one or both of the highly conserved prolyl 
residues located in ODDD of HIF-a (HIF-la and HIF- 
2a) becomes hydroxylated by the prolyl-4-hydroxylase 
(PHD), and then the hydroxylated HIF-a interacts with 
the von-Hippel-Lindau tumor suppressor (VHL) and re- 
cruits ubiquitin ligase, which is consequentiy degraded by 
the proteasome [19,22,30,31]. However, under hypoxic 
conditions, HIF-a, which is stabilized, translocates to the 
nucleus via its nuclear localization signal (NLS) motif, 
where it dimerizes with HIF-|3 into HIF heterodimer and 
binds to the core DNA motif (G/ACGTG) in hypoxia- 
response elements (HREs) for transcriptional activation of 
the target genes [19,22,32,33]. 

The teleost-specific whole-genome duplication that oc- 
curred early in the teleost evolution generated six hif-a 
genes: 1A/B, 2A/B, and 3A/B, which was followed by the 
loss of one member of each A/B paralogous pair in most 
euteleosts but not in cyprinids [20]. The schizothoracine 
fish (Teleostei: Cyprinidae), which are divided into three 
groups (primitive, specialized, and highly specialized) ac- 
cording to the differences in their scales, pharyngeal 
teeth, and barbels, are the endemic and most diverse 
group of cyprinids distributed in the Tibetan Plateau and 
its adjacent areas. The elevation in this area ranges from 
700 to 5000 m, and the environment is characterized by 
hypoxia and low temperatures [34-36]. Primitive schi- 
zothoracine fish including Schizothorax, Schizocypris 
(not found in China), and Aspiorhynchus, are distributed 
at the edge of the Tibetan Plateau. Specialized schi- 
zothoracine fish included Diptychus, Ptychobarbus, and 
Gymnodiptychus, and highly specialized schizothoracine 
fish included Gymnocypris, Oxygymnocypris, Schizopy- 
gopsis, Platypharodon, Chuanchia, and Herzensteinia. 
The distribution of specialized and highly specialized 
schizothoracine fish is largely limited to the central re- 
gion of the plateau [37]. Each group represents specific 
periods of geological evolution of the Tibetan Plateau, 



which ultimately caused the modern distribution pat- 
tern of schizothoracine fish, in which primitive, special- 
ized, and highly specialized schizothoracine fish are 
distribued from the edges to the central part of the 
plateau (Figure 1) [35]. Previous analyses of the mito- 
chondrial genome have demonstrated that schizothora- 
cine fish, which are teleost fishes, are well-adapted to 
high altitudes and hypoxia [38]. However, hif-a genes, 
which are related to hypoxia, have been little studied in 
schizothoracine fish. Chi [12] referred to the evolution- 
ary patterns of the hif-a genes in two high-altitude fish, 
but no signs of adaptive evolution was detected since 
the species limitation. 

To better understand the relationship between adapta- 
tion to hypoxia and hif-a genes in plateau and plains 
fish, teleost-specific hif-a gene duplications (hif-laAlB 
and hif-2aAIB) were cloned and identified in six species 
of schizothoracinae and in zebrafish. Further computa- 
tional and experimental analyses reveal the mode of evo- 
lution and response to hypoxia of these gene paralogs. 

Results 

Characterization of hif-a duplications and phylogenetic 
analysis 

Although many studies have focused on HIF in cyprinids, 
duplicated HIF-a paralogs have never been identified in 
schizothoracine fish. The sequences of hif-a paralogs {hif- 
laA/B and hif-2aA/B) of the six species of schizothoraci- 
nae were obtained here for the first time [GenBank: 
KJ679876-KJ679899]. The N-terminal of zebrafish hif-laA 
sequence was also isolated [GenBank: KJ679875]. The data 
sets supporting the results of this article are also available 
in the Dryad Digital repository [39]. The lengths of differ- 
ent hif-a paralogs in schizothoracine fish were found to 
vary. However, the lengths of specific hif-a genes remained 
consistent among highly specialized schizothoracine fish, 
even when the schizothoracine fish came from different 
genera (Additional file 1: Table SI). 

Protein analysis showed that HIF-a paralogs from 
schizothoracine fish all contained the main domains, like 
their mammalian counterparts. The similarity between 
duplicated HIF-laA and HIF-laB was 52%, and the 
similarity between duplicated HIF-2aA and HIF-2aB 
was 55%. The bHLH-PAS domain of HIF-as N-terminal 
was highly conserved, but the C-terminal was less con- 
served, especially near the oxygen-dependent proline hy- 
droxylation sites. A deletion, specifically, the conserved 
proline hydroxylation motif LxxLAP, was detected in the 
NODD domain of HIF-laA in all six schizothoracine fish 
and zebrafish (Figure 2 and Additional file 2: Figure SI A). 
Currently, hif-laA was only amplified in Cyprinidae, so 
this deletion might be Cyprinidae-specific. Sequence align- 
ment suggested that there was a schizothoracine-specific 
mutation in the CODD domain of the HIF-laB of 
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Figure 1 Altitude distribution of three groups of schizothoracine fish in Tibetan Plateau. The map was downloaded from the National 

Geographic ArcMap Esri db on June 22, 2014 [79]. The solid boxes indicate the altitude distribution range of three groups of the main Schizothoracine 

genera. I: primitive schizothoracine fish, II: specialized schizothoracine fish, III: highly specialized schizothoracine fish. 
* ) 



specialized and highly specialized schizothoracine fish, 
specifically the proline hydroxylation motif mutated into 
PxxLAP (Figure 2 and Additional file 2: Figure SI B). All 
deletions and mutations were located in functional ODD 
domains that interacted with PHD via the PHD-HIF 
oxygen-sensing system [11,22,40]. This may affect func- 
tion during biological processes. The targets of PHD in 
mammals are proline 402 and proline 564 in the ODD do- 
main, but Pro-564 is the primary critical substrate in 
PHD binding to HIF-la [10,41]. It is here speculated 
that HIF-locB is more likely to change function than 
other HIF-cts are. 

The chordate phylogenetic tree of hif-la and hif-2a 
genes was reconstructed using the MrBayes and RAxML 



softwares with sea squirt and amphioxus hif-a sequences 
as the outgroups [42,43]. The two computer programs 
inferred similar tree topologies. The Bayesian posterior 
probability values and maximum likelihood bootstrap 
values were high (Figure 3, Additional file 3: Figure S2). 
The topology was largely consistent with that found in 
previous studies [11,29]. However, the phylogeny esti- 
mated for the hif-la gene was inconsistent with the 
recently published hif-la phylogenies based on the Phy- 
loBayes [44] and RAxML softwares [20]. Recent results 
indicated a sister relationship between hif-laA and 
hif-laB, our results indicate instead that hif-laB might 
be more similar to the ancient hif-la than hif-laA is. 
Massey ef al. [45] and Roje [46] have researched on the 
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Figure 2 Partial multiple sequence alignment of the deduced HIF-1a protein sequences. A) In Cyprinidae, HIF-laA proteins were found to 
have a deletion (LxxLAP) in the conserved NODD domain. B) Specialized and highly specialized schizothoracine fish were found to have a 
mutation (LxxLAP to PxxLAP) in the conserved CODD domain of HIF-laB. Dashes indicate that the gaps inserted for facilitated alignment. The 
two conserved proline hydroxylation motif LxxLAP are indicated by the dotted boxes. 



effects of non-neutral gene on phylogeny. In this study, 
more hif-laB sequences were used than any other stud- 
ies, which might have led to a non-synonymous to 



synonymous rate ratio (dN/dS) distinct from that ob- 
tained in previous studies. The phylogeny (Additional 
file 4: Figure S3 A) in the schizothoracine fish constructed 
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Figure 3 Phylogenetic tree was constructed for the selected chordate hif-la and hif-2a genes. The phylogenetic tree was constructed 
using MrBayes with GTR + 1 (0.0821) + G (1.0813) modes. To search tree space, Markov chain Monte Carlo (MCMC) simulations were performed for 
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MCMC analysis. The Bayesian posterior probability values are indicated beside the branches. 
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using concatenated hif-la and hif-2a sequences agrees 
with that using mitochondrial cytochrome b [38,47]. For 
this reason, it was considered suitable for use in the fol- 
lowing selection pressures analysis. 

Molecular evolution of the teleostei HIF-a duplicate 

First, standard branch-site models in PAML [48] were 
used to assess the different selective pressures experi- 
enced by the protein coding sequences of each HIF-a 
paralogous pair. In this model, the foreground branches 
are the branches tested for positive selection, and the 
background branches are all other branches on the tree 
[48]. In the likelihood ratio test (LRT), branch-site model 
A is the alternative model allowing CO2 to vary among co- 
dons, while the simpler null model is model A but with 
(1)2 = 1 fixed [48]. In this study, three methods were used 
to specify the foreground branches. According to the 
LRT of branch-site models, for HIF-laB, when the clade 
of all schizothoracine fish was specified as the foreground 
branches, the likelihood obtained under model A was 
not significantly higher than that obtained under the null 
hypothesis of no positive selection; but when the clade 
of specialized and highly specialized schizothoracine fish 
was specified as the foreground branches, the likelihood 
obtained under model A was significantly higher than 
that of the null model (2AlnL = 12.16, p < 0.01, df= 1) 
and go > 1 (w = 7.11, Table 1). When only the clade of 
highly specialized schizothoracine fish was specified as 
the foreground branches, the likelihood obtained under 
model A was quite significantly higher than that of the 
null model (2AlnL = 19.79, p < 0.01, df= 1) and 0 > 1 
(co = 23.24) (Table 1). The posterior probability of sites 
that had undergone positive selection was calculated 
using the Bayes empirical Bayes (BEB) method in 
model A. For these last two analyses, there were 9 and 
10 amino acid sites in HIF-laB sequences with a pos- 
terior probability greater than 0.5. In the second 
method, the amino acid sites at position 464 (accord- 
ing to HIF-laB of D. rerio) had a posterior probability 
0.990 as indicated by BEB, which was significant at the 
1% level, but in the third method, the amino acid sites 
at position 387 and 510 had a posterior probability 
0.989 and 0.987 as indicated by BEB, which was signifi- 
cant at the 5% level. In this way, they were identified 
as crucial amino acid sites that had experienced 



positive selection. Three of the 9 sites and 4 of the 10 
sites were located in the ODD domain by the last two 
methods. For HIF-laA, HIF-2aA, and HIF-2aB, no 
other positively selected site was detected. 

Secondly, stochastic branch-site models in fitmodel 
[49] were used to detect the different selective pressures 
experienced on the protein coding sequences of HIF-la 
and HIF-2a, respectively. For HIF-la, the likelihood ob- 
tained under M3 + SI (co2 = 8.7857) model was signifi- 
cantly higher than that obtained under null (Ml + SI) 
model with co2 > 1 and 2AlnL = 6.29 (p < 0.05). However, 
there were only 3 amino acid sites (77 K, 395S, and 
398 L, according to HIF-laB of D. rerio) with a posterior 
probability > 0.5. For HIF-2a, no other positively selected 
site was detected. 

In standard branch-site models, although the fore- 
ground branches specified for HIF-laB were different, 
the last two methods showed similar positively selected 
sites with different levels of significance. With stochastic 
branch-site models, positive selection was also detected 
in HIF-laA and HIF-laB. It is here hypothesized that 
the HIF-laB experienced strong selective pressure, 
which may have led to changes in function. The higher 
the altitude at which the schizothoracine fish were dis- 
tributed, the more significant the selective pressure was. 
This may be one of the mechanisms by which schi- 
zothoracine fish have adapted to the high altitudes of the 
Tibetan Plateau. 

Expression and hypoxic regulation of different HIF-a 
paralogs in HEK 293 T cells 

RytKonen [20] showed that the transcription of hif-la de- 
creased after hypoxic insult in adult zebrafish, but the 
transcription of hif-2aA remained stable and hif-2aB tran- 
scription increased significantly. To determine whether 
changes in oxygen tension regulate HIF-a at the protein 
level, HEK 293 T cells were transfected with either the in- 
dicated Myc-tagged HIF-a constructs or an empty vector, 
and cells were subjected to 18 h of hypoxic treatment after 
6 h of normoxic cultivation. Western blot analysis indi- 
cated that all the fish HIF-as could be expressed in HEK 
293 T cells. Compared to the normoxic condition, there 
were only rarely significant increases in levels of protein 
expression under hypoxic conditions: HIF-laA and 
HIF-laB in D. rerio, HIF-laB in Gd. pachycheilus and 



Table 1 Positively selected sites detected in HIF-laB of specialized and highly specialized schizothoracine fish by 
branch-site models 



Foreground branch 


Aln L Parameter estimates 


Positive sites 


specialized and highly specialized 


6.08 pa = 0.820 p, = 01 49 p 2 = 0.026 0J 0 = 0.074 oj, = 1 oj 2 = 7.1 1 3 


21 5P 370E 371 T 386P 387E 396S 464D** 


schizothoracine fish 




510D 534I 


highly specialized schizothoracine fish 


9.89 p 0 = 0.820 p, = 01 7 p 2 = 0.01 0 a> 0 = 0.079 w, = 1 oj 2 = 23.246 


41 H 370E 371 T 386P 387E* 428S 464D 






467P 510D* 648 T 



Note: * mean pp > 95%, ** mean pp > 99%. 
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HIF-2aA of G. eckloni (Figure 4). Generally, the expression 
level of HIF-loc was greater than that of HIF-2oc after hyp- 
oxic treatment and under normoxic conditions (Additional 
file 5: Figure S4). 

To determine whether selective pressure could have 
affected the transcriptional activity of HIF-a toward their 
common targets, empty vector (control) and the indi- 
cated Myc-tagged HIF-a constructs were co-transfected 
with the HRE-luciferase reporter and pTK-Renilla lucif- 
erase in HEK 293 T cells. Each HIF-la and HIF-2a was 
found to upregulate luciferase activity under normoxic 
conditions (Figure 5). This indicated that the function of 
fish HIF-as had been conserved. In addition, HIF-laB 
showed more luciferase activity than other HIF-as, even 
under normoxic conditions, suggesting that HIF-laB 
protein may be more stable. The transcriptional activity 
of HIF-2a was obvious under normoxic conditions but 
not detectable under hypoxic conditions. 

Discussion 

Fish possess homologs of HIF-a and HIF- [3 that share 
strong similarity with those of humans. They may play 
roles similar to those of their mammalian counterparts 
in expression of oxgen-dependent genes [13]. The se- 
quence of hif-a gene had been extensively characterized 
in fishes. The first fish hif-a gene was identified in 



rainbow trout {Oncorhynchus mykiss) in 2001 [50]. The 
sequence was slightly shorter than that of mammals. It 
included bHLH-PAS and ODD domains and proline and 
asparagine residues, which were relatively conserved in 
the deduced rainbow trout HIF-a protein. Subsequently, 
hif-a gene was reported from killfish (Fundulus hetero- 
clitus), grass carp {Ctenopharyngodon idellus), crucian 
carp (Carassius carassius), zebrafish (Danio rerio), Atlantic 
croaker (Micropogonias undulatus), asp {Aspius aspius), 
three-spined stickleback (Gasterosteus aculeatus), Russian 
sturgeon {Acipencer gueldenstaedtii), naked carp {Gym- 
nocypris przewalskii), Wuchang bream (Megalobrama 
amblycephala), Chinese sucker (Myxocyprinus asiaticus), 
bighead carp {Hypophthalmichthys nobilis), silver carp 
(Hypophthalmichthys molitrix), Prenant's schizothoracin 
(Schizothomx prenanti), Namucuo naked carp (Gymnocy- 
pris namensis), and so on [10,12,51-58]. 

Although previous studies identified HIF-a in schi- 
zothoracine fish, the present study is the first to identify 
teleost-specific duplicated HIF-a paralogs. It not only 
provides new information for HIF research in teleosts 
but also supplies a foundation for further study of the 
adaptation of schizothoracine fish to high altitudes and 
low oxygen levels in the Tibetan Plateau. Each schi- 
zothoracine fish HIF-a was smaller than the correspond- 
ing human HIF-a. The HIF-laA/B and HIF-2aA/B of 
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Figure 4 Expression of HIF-a in HEK 293 T cells under normoxic and hypoxic conditions. N: normoxia, H: hypoxia. There was no significant 
different in protein expression levels under hypoxic and normoxic conditons, excepting that there was more expression of HIF-1 aA and HIF-1 aB 
of D. rerio, HIF-laB of Gd. Pachycheilus and HIF-2aA of G. eckloni under hypoxic conditions. 
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Figure 5 Transcriptional activity of HIF-as under normoxic and hypoxic conditions. Each HIF-1 a and HIF-2ct was found to upregulate 
luciferase activity under normoxia, and the luciferase activity of HIF-1 aB was high even under normoxia. 



schizothoracine fish shared more than 46 and 53% 
similarity with human HIF-loc and HIF-2a, respectively. 
HIF-laB and HIF-2aB shared more sequence identities 
with human HIF-lct and HIF-2oc than HIF-laA and 
HIF-2aA did in all species except Gd. pachycheilus. 
The bHLH-PAS domain of HIF-as N-terminal was 
highly conserved, but the ODD domain of the C- 
terminal was less conserved, especially in regions in the 
vicinity of the oxygen-dependent proline hydroxylation 
sites. A Cyprinidae-specific deletion (LxxLAP) in the 
NODD domain of HIF-laA may have occurred in all 
six schizothoracine fish, and a schizothoracine-specific 
mutation LxxLAP mutated to PxxLAP in the CODD 
domain of the specialized and highly specialized schi- 
zothoracine fish HIF-laB. The HIF hydroxylase pathway 
was regulated by the PHD family of oxygen-dependent 
prolyl hydroxylases in metazoans [19]. In mammals, the 
Pro-564 in the CODD domain is the primary critical sub- 
strate in PHD binding to HIF-la [31]. It was here deduced 
that although the deletion and mutation were all located 
in the ODD domain, the change in the CODD domain of 
HIF-laB was more important than the deletion in the 
NODD domain of HIF-laA was in the evolution of schi- 
zothoracine fish. 

The HIF pathway plays a pivotal role in the response 
to hypoxia, and the PHD family of oxygen-dependent 



prolyl hydroxylases plays a critical role in regulating HIF 
stability [19]. Some candidate genes associated with the 
HIF pathway showed signals of positive selection in spe- 
cies living in the Tibetan Plateau. For instance, EPAS1 
(HIF-2a), EGLN1, and PPARA were detected in Tibetans 
[59,60]; ADAM17 in the yak [61]; ADORA2A, CCL2, 
ENG, PIK3C2A, PKLR, ATP12A, and NOS3 in the Ti- 
betan antelope [62]; SRF, TXNRD2, and WNT7B in the 
ground tit [63]; and EPAS1, SIRT7, PLXNA4, and MAFG 
in the Tibetan mastiff [64]. In this study, computational 
estimation suggested that the specialized and highly spe- 
cialized schizothoracine fish HIF-laB have experienced 
significantly selective pressure. Positively selected sites 
were detected in HIF-la, HIF-2a, and HIF-2a, respectively 
in schizothoracine fish, Tibetans, and Tibetan mastiff. 
This shows that fish and mammals developed different 
mechanisms for adapting to the special environment of 
the Tibetan Plateau. 

Functional analysis was performed by transfecting 
schizothoracine fish HIF-a in HEK 293 T cells. Results 
showed that the expression level of HIF-la differed 
from that of HIF-2a. Hypoxia was found to increase 
the abundance of HIF-laA/B but not that of HIF- 
2aA/B. The decrease in the protein levels of HIF-2aA/B 
may be attributable to the combined impact of the re- 
duced protein synthesis under hypoxic conditions and 
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the rapid folding rate of HIF-oc protein in mammalian 
cells [65]. The luciferase assay showed that all the HIF-ocs 
characterized in this study were capable of forming func- 
tional heterodimers with human HIF-(3 and could acti- 
vate HRE reporter gene, but they performed different 
levels of oxygen-dependent regulation. Each HIF-a (HIF- 
laA/B and HIF-2aA/B) was able to upregulate luciferase 
activity under normoxic conditions, and HIF-laB was 
more significant, suggesting that HIF-laB was more 
stable under normoxic conditions. The transcriptional 
activity of HIF-laB was higher than that of other HIF- 
las under hypoxic conditions, and transcriptional activity 
was stronger under hypoxic conditions than under nor- 
moxic conditions. These results were also supported by 
structural analysis. For HIF-laB, the bHLH-PAS domain, 
which is responsible for heterodimerization with HIF-(3 
and DNA binding, was conserved; and the ODD domain, 
which interacts with PHD, was mutated [10,19,22,66,67]. 
It was also discovered that HIF-2aA/B downregulated 
luciferase activity under hypoxic conditions. This was 
consistent with the decrease observed in protein levels 
of HIF-2aA/B after hypoxia treatment. 

Conclusions 

The teleost-specific duplicated hif-a paralogs, hif-laA, 
hif-laB, hif-2aA, and hif-2aB, were successfully identi- 
fied in schizothoracine fish. Each deduced HIF-a was 
found to have the same principal domains as its mam- 
malian counterparts, and had evolved specialized roles 
in the response to hypoxia. It is here postulated that 
HIF-laB may be the most important regulator in the 
adaptation of schizothoracine fish to the environment of 
the Tibetan Plateau. Although sequence variation and 
significant selective pressure were detected in HIF-laB 
of the specialized and highly specialized schizothoracine 
fish, transcriptional activity showed no significant differ- 
ence between plateau fish (schizothoracine fish) and 
plains fish (zebrafish). This might be because zebrafish 
are also hypoxia tolerant fish [10,13]. Further studies are 
needed to explore the structure and functions of schi- 
zothoracine fish HIF-as and to determine their roles in 
the mechanism underlying the response to hypoxic condi- 
tions. The adaptation to Tibetan Plateau of schizothora- 
cine fish involves a complicated evolutionary mechanism, 
which is affected by a cluster of genes. In the future, re- 
search on transcriptome and genome of schizothoracine 
fish may contribute to explain the mechanism of their 
adaption to Tibetan Plateau. 

Methods 

Ethics statement 

The experiments were performed in accordance with the 
Ethics Committee of the Institute of Hydrobiology, 
Chinese Academy of Sciences. The policies were enacted 



according to Chinese Association for Laboratory Animal 
Sciences, and coordinated with the Institutional Animal 
Care and Use Committee (IACUC) protocols [68,69]. 

Sampling 

Schizothorax prenanti (primitive schizothoracine fish) and 
Gymnodiptychus pachycheilus (specialized schizothoracine 
fish) were collected from Ya'an in Sichuan Province; Gym- 
nocypris eckloni, Schizopygopsis pylzovi, and Platypharo- 
don extremus (highly specialized schizothoracine fish) 
were collected from the headwaters of the Yellow River; 
Gymnocypris namensis (highly specialized schizothoracine 
fish) was collected from Namtso Lake; Danio rerio (AB 
strain) was collected from Institute of Hydrobiology, 
Chinese Academy of Sciences. Permissions for extract- 
ing the fish were obtained from Aquaculture Bureau of 
Sichuan Province and Aquaculture Bureau of Qinghai 
Province. The seven species, which are all the hypoxia- 
tolerant fish, live at different altitudes. The schizothora- 
cine fish cover five genera and three groups. Tissue sam- 
ples (heart and liver) were obtained and immediately 
infused with Trizol (TaKaRa) and conserved at -70°C until 
RNA extraction. 

RNA isolation, cDNA synthesis, and sequence 
amplification 

RNA extraction was performed using the Trizol reagent 
according to the manufacturer's instructions. The quality 
of RNA was determined on 1.2% EtBr- agarose gels, and 
the quantity (A 2 6o/A 2 so) was measured with a UV spec- 
trophotometer (Thermo Scientific Nanodrop 2000). 

First, 1 ug of total RNA was used to synthesize first- 
strand cDNA using reverse transcriptase (Promega) in 
the presence of RNAse inhibitor (Invitrogen) with a 
cocktail of oligoT and dNTP (TaKaRa). The 5 -RACE- 
ready cDNA and 3 -RACE-ready cDNA were synthe- 
sized using a SMART™ RACE cDNA Amplification Kit 
(ClonTech) and FirstChoice® RLM-RACE Kit (Ambion), 
respectively. To obtain the intact cDNA, different tissue 
samples with the best quality RNA were chosen: heart 
for Schizothorax prenanti, Gymnodiptychus pachychei- 
lus, and Gymnocypris eckloni, and liver for Gymnocypris 
namensis, Schizopygopsis pylzovi, Platypharodon extre- 
mus, and Danio rerio. 

Universal primers (Additional file 6: Table S2) for ampli- 
fying partial cDNA fragments were designed using the 
Primer Premier 5 program against the conserved consen- 
sus sequences derived from multiple alignment of 
hif-laA/B and hif-2aA/B sequences from Schizothorax 
prenanti (Prenant's schizothoracin), Danio rerio (zebra- 
fish), Hypophthalmichthys molitrix (silver carp), and 
Hypophthalmichthys nobilis (bighead carp), which were 
collected from Ensembl or National Center for Biotech- 
nology Information (NCBI, Additional file 7: Table S3) [70]. 
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First-strand cDNAs were used in the primary PCR 
amplification with the universal primers. The condi- 
tions for PCR amplification were an initial denaturation 
at 94°C for 3 min, 25 cycles of denaturation at 94°C for 
30 s, annealing at 60°C for 30 s, and extension at 72°C 
for 2 min, and a final extension at 72°C for 10 min. 

The partial sequences obtained with the universal 
primers were used to design the gene specific primers for 
5'- and 3 '-RACE reactions (Additional file 8: Table S4). 
RACE reactions were performed in accordance with the 
instructions in a SMART'" RACE cDNA Amplification Kit 
(ClonTech) and FirstChoice 8 RLM-RACE Kit (Ambion). 

All PCR products were detected in 0.8% EtBr-agarose gels, 
and bands corresponding to hif fragments were excised 
and extracted with a DNA Gel Extraction Kit (Tiangen). 
The sequences were recombined with pMD18-T vector 
(TaKaRa), sequenced using ABI3730XL sequencer (Beijing 
Tianyi Huiyuan Bioscience and Technology Inc.), and 
assembled using SeqMan in DNASTAR Lasergene v7.1 
software (DNASTAR Inc.). The assembled sequences 
were verified using the NCBI database. The full-length 
sequences of hif-laA/B and hif-2aA/B were amplified 
using newly designed gene-specific primers to the 5' 
and 3' ends with the restriction enzyme sites (Additional 
file 9: Table S5). 

Sequence alignment and phylogenetic analysis 

The alignment of the hif-la and hif-2a coding sequences 
(corresponding to a total of 69 sequences and 3,090 bases 
pairs obtained after alignment) was based on notional 
translations of nucleotide sequences using ClustalX [71]. 
Additionally, the GTR + I (0.0821) + G (1.0813) model was 
adopted for the phylogeny by running PAUP4.0bl0 and 
modeltest3.7 [72,73]. Finally, Mrbayes-3.1.2 with 2,000,000 
iterations and RAxML with 1,000 nonparametric boot- 
strap replicates were used to construct a phylogenetic tree 
[42,43]. The sequences of sea squirt and amphioxus hif-a 
were used as outgroups. 

Model testing of selection pressures 

To determine whether hif gene paralogs have under- 
gone statistically significant differences in selection 
pressures, standard branch-site models in codeml soft- 
ware of PAML (version 4.6) and stochastic branch-site 
models in fitmodel (version 20140407) were used 
[48,49,74-78]. The two methods all employed ML esti- 
mates of the ratio of nonsynonymous to synonymous 
substitutions (dfj/ds = to) and nested likelihood ratio tests 
(LRTs) on a phylogeny tree (Additional file 4: Figure S3). 
In standard branch-site models, only the Cyprinidae spe- 
cies which had all the four paralogs were included, and 
the clades corresponding to schizothoracine fish were 
marked as foreground branches. 



Recombined plasmid construction 

The correct and intact hif-la and hif-2a genes of Danio 
rerio, Schizothorax prenanti, Gymnodiptychus pachy- 
cheilus, and Gymnocypris eckloni were isolated from 
one of the TA clones and subcloned into pCMV-Myc 
vector (CloneTech) using special restriction enzyme 
sites (Additional file 6: Table S2). All recombined plas- 
mids were verified by DNA sequencing. 

Cell culture 

Human embryonic kidney (HEK) 293 T cells were cul- 
tured in Dulbecco's modified Eagle's medium (DMEM, 
Hyclone) supplemented with 10% fetal bovine serum 
(FBS, Hyclone) in a 5% CO2 incubator (Thermo). Hyp- 
oxia treatment was conducted in a hypoxia carbon diox- 
ide incubator filled with 1% O2. 

Western blot analysis 

HEK 293 T cells were seeded in 6-well plates for 24 h 
before transfection, and then transfected with 1.6 ug 
every plasmid (included empty vector and indicated 
Myc-tagged HIF-a constructs) using VigoFect (Vigor- 
ous). Six hours after transfection, the cells were divided 
into two groups and cultured under normoxic and hyp- 
oxic conditions (1% O2) for 18 h. Cells were lysed in 
RIPA buffer (Beyotime) and 1% PMSF was added (Beyo- 
time). Then the lysate was mixed with 5x protein load- 
ing buffer (Beyotime), and boiled for 5 min. Protein 
samples were separated using the denaturing electro- 
phoresis on 12% SDS-polyacrylamide vertical gels and 
transferred to a PVDF (polyvinylidene fluoride) mem- 
brane (Millipore). The membrane was blocked for 2 h in 
5% nonfat dry milk in TBST (TBS: 7.3 g NaCl and 
3.03 g Tris-Base diluted in 1,000 mL H z O, TBST: TBS 
supplemented with 0.1% Tween 20) at room 
temperature and incubated with primary antibody (c- 
Myc (9E10), Santa Cruz) that had been diluted 1:3,000 
in 1% nonfat dry milk in TBST for 12 h at 4°C. The 
membrane was washed three times with TBST for 8 min 
each and incubated with the secondary antibody (HRP- 
labeled Goat Anti-Mouse IgG (H + L), Beyotime) that 
had been diluted 1:2,000 in 1% nonfat milk in TBST. 
The membrane was then washed, and the HIF-a protein 
was detected using enhanced chemiluminescence (ECL, 
Millipore) in an ImageQuant™ LAS 4000 mini (GE Health- 
care). Tubulin was used as the internal control. 

Luciferase detection 

HEK 293 T cells were seeded in 24-well plates for 24 h be- 
fore transfection and then co-transfected with empty vector 
(control) or the indicated Myc-tagged HIF-a constructs, 
the HRE-luciferase reporter (provided by Professor Wuhan 
Xiao), and the pTK-Renilla luciferase reporter (used as an 
internal control) using VigoFect (Vigorous). Six hours after 
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transfection, the cells were separated into two groups and 
cultured under normoxic and hypoxic conditions (1% O2) 
for 18 h. Luciferase activity was measured using a Dual- 
Luciferase Reporter Assay System (Promega). 

Additional files 
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